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LONG-TERM  GOALS 

The  long-term  goal  of  this  research  is  to  develop  optimal  array  signal  processing  techniques  for  passive 
source  localization  in  littoral  shallow-water  environments.  It  has  long  been  known  that  multi-modal 
dispersion  in  a  shallow  water  waveguide  degrades  the  performance  of  bearing  estimates  by 
conventional  plane -wave  beamforming.  This  is  due  to  spurious  effects  unique  to  the  waveguide 
environment  such  as  multiple  peaks  and  beam  spreading  in  the  beamformer  output  [1,2].  We  have 
developed  array  signal  processing  techniques  that  account  for  and  exploit  multi-modal  wave 
propagation  and  dispersion. 

OBJECTIVES 

The  primary  objectives  of  the  research  are  to: 

■  Develop  the  array  invariant  passive  localization  techniques  that  require  little  a  priori 
knowledge  of  the  wave  propagation  environment. 

■  Instantaneously  and  simultaneously  localize  multiple  sources  in  a  shallow-water  waveguide 
without  extensive  computations  and  without  ambiguity. 

APPROACH 

The  array  invariant  method  for  passive  localization  of  impulsive  sources  has  been  theoretically  derived 
and  experimentally  demonstrated  in  FY05  effort  using  the  data  acquired  fonn  the  Acoustic  Clutter 
Experiment  in  2001  and  2003.  It  has  been  shown  that  simple  and  robust  passive  source  localization  is 
possible  using  the  array  invariant  method  without  a  priori  knowledge  of  environmental  parameters.  We 
have  extended  the  array  invariant  method  for  range  and  bearing  estimation  of  a  broadband  random 
noise  source  that  is  not  necessarily  impulsive  in  the  time  domain,  while  maintaining  the  advantages  of 
the  array  invariant  method  for  impulsive  source  localization.  We  have  theoretically  and  numerically 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2007 

2.  REPORT  TYPE 

Annual 

3.  DATES  COVERED 

00-00-2007  to  00-00-2007 

4.  TITLE  AND  SUBTITLE 

ONR  Graduate  Traineeship  Award  In  Ocean  Acoustics  For  Sunwoong 

Lee 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology, 77  Massachusetts  Avenue,  Room 
5-212, Cambridge, MA, 02139 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

code  1  only 

14.  ABSTRACT 

The  long-term  goal  of  this  research  is  to  develop  optimal  array  signal  processing  techniques  for  passive 
source  localization  in  littoral  shallow-water  environments.  It  has  long  been  known  that  multi-modal 
dispersion  in  a  shallow  water  waveguide  degrades  the  performance  of  bearing  estimates  by  conventional 
plane-wave  beamforming.  This  is  due  to  spurious  effects  unique  to  the  waveguide  environment  such  as 
multiple  peaks  and  beam  spreading  in  the  beamformer  output  [1,  2].  We  have  developed  array  signal 
processing  techniques  that  account  for  and  exploit  multi-modal  wave  propagation  and  dispersion. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


shown  that  the  array  invariant  method  can  be  used  not  only  for  localization  of  a  single  source,  but  for 
simultaneous  localization  of  multiple  noise  sources  in  an  ocean  waveguide  without  ambiguity. 

WORK  COMPLETED 

The  research  work  in  Fiscal  Year  2006  was  a  great  success  in  extending  the  array  invariant  source 
localization  techniques  in  a  dispersive  waveguide  that  require  little  a  priori  knowledge  of  the 
environment.  The  extended  array  invariant  techniques  enable  instantaneous  and  simultaneous 
localization  of  multiple  random  noise  sources  in  a  horizontally-stratified  ocean  waveguide  from 
passive  beam-time  intensity  data  obtained  after  conventional  plane-wave  beamforming  of  acoustic 
array  measurements. 

Localization  of  multiple  sources  using  the  array  invariant  method  has  advantages  over  the  conventional 
triangulation  method  using  two  arrays.  First,  the  array  invariant  method  does  not  require  two  arrays 
with  sufficient  spatial  separation,  which  is  a  significant  advantage  when  the  localization  is  made  from  a 
mobile  platform.  Second,  the  array  invariant  method  does  not  suffer  from  the  ghost  source  problem 
typical  in  the  triangulation  method  when  there  are  multiple  sources  present  [3].  The  array  invariant 
method  is  significantly  advantageous  over  MFP,  where  unambiguous  source  range  estimation  is  nearly 
impossible  even  with  accurate  environmental  knowledge,  when  there  are  unknown  number  of  multiple, 
spatially  stationary  sources  in  an  ocean  waveguide  [4,  5]. 

RESULTS 

We  have  previously  shown  that  the  localization  of  an  impulsive  source  can  be  achieved  using  the  array 
invariant  method  without  requiring  the  a  priori  knowledge  of  environmental  parameters  and  without 
extensive  computations  [6].  We  have  further  extended  the  method  to  show  that  the  array  invariant 
method  can  be  applied  for  localization  of  a  broadband  random  noise  sources.  This  is  achieved  by 
cross-correlating  the  beam-time  series  with  the  time  series  measured  by  the  acoustic  sensor  at  the 
center  of  the  array  [7]. 

We  first  show  that  the  range  and  bearing  of  a  broadband  random  noise  source  can  be  determined  using 
the  extended  array  invariant  method.  The  source  is  located  at  ro  =5  km  and  9o  =  60°  as  shown  in  Fig. 
1.1.  The  environment  is  100-m  deep  Pekeris  waveguide  with  sand  bottom.  The  source  and  the 
receiving  array  depths  are  50  m  and  30  m,  respectively.  The  aperture  of  the  receiving  array  is  150  m. 
The  source  power  spectral  density  is  assumed  to  be  0  dB  re  1  pPa2  /  Hz  within  the  390  to  440-Hz 
frequency  band. 
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Figure  1:  Top  view  of  the  geometry  for  the  single  source  localization  example.  The  source  is  at  60° 
from  the  broadside  of  the  horizontal  receiving  array,  and  at  5-km  range. 


(a)  Cross-correlation  (b)  Beamformer  migration 


Scan  Angle  (s  =  sin  0)  Scan  Angle  (s  =  sin  0) 

Figure  2:  The  cross-correlated  intensity  field  for  the  source-receiver  geometry  shown  in  Fig.  1  in  the 
sand-bottom  Pekeris  waveguide  environment.  The  black  solid  line  overlain  in  the  r  >  0  domain  is  the 
beamformer  migration  line  in  terms  of  the  reduced  travel  time.  The  black  solid  line  in  the  r  <0 

domain  is  at  sind . 


The  cross-correlated  intensity  field  IBo  (s,x  )  is  shown  in  Fig.  2,  where  s  =  sin9  is  the  array  scan  angle, 
and  x  is  the  delay  time.  Two  black  solid  lines  are  overlain  on  Fig.  2.  One  line  in  the  x  >  0  domain  is  the 
beamformer  migration  line  ~  s  (x  )  for  an  impulsive  source  in  an  ideal  waveguide  as  defined  in  Eq.  (8) 
of  Ref.  [6],  and  the  other  line  in  the  x  <  0  domain  marks  the  source  bearing  sinOo  .  These  two  lines 
show  that  the  cross  correlated  intensity  field  is  bounded  by  these  lines  except  that  the  field  is  smeared 
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by  the  beampattern  of  the  horizontal  receiving  line  array.  Given  these  two  bounding  lines  in  the  cross- 
correlated  intensity  field,  we  can  estimate  both  source  range  and  bearing  without  ambiguity. 

In  order  to  identify  these  two  bounding  lines  from  Fig.  2,  we  apply  an  image  transform  technique 
similar  to  the  Radon  transform  [8]  to  the  cross-correlated  beam-time  image.  In  this  transform,  a  given 
image  is  integrated  along  a  semi-infinite  line  starting  from  (s,  0)  and  having  an  angle  (pr  from  the 

positive  x  -axis  in  the  clockwise  direction.  Now  we  show  that  both  the  range  and  bearing  of  the  source 
can  be  extracted  from  the  transformed  intensity  field  in  Fig.  3. 
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Figure  3:  The  transformed  intensity  image  In  (s,(pr  ) 
of  the  cross-correlated  intensity  image  in  Fig.  2. 

First,  the  bearing  of  the  source  can  be  estimated  from  the  value  of  $s$  where  the  maximum  of  In  (s,cpr 
=  180°)  occurs.  In  Fig.  3,  this  peak  is  seen  to  occur  at  s/'o=  0.853  or  0"  o  =  58.5°  .  This  is  within  3\% 
of  the  true  source  bearing  0o  =  60°  .  Second,  once  the  bearing  of  the  source  has  been  estimated,  the 
range  of  the  source  can  be  estimated  using  c(z)sin9''  o  r'o  =—  ,  tamp"  r  where  cp"  r  =  arg  maxcp  In  (s  = 
s^cpr  ).  In  Fig.  3,  this  peak  occurs  at  qF  r=  346.7°  .  The  range  of  the  source  is  r  then  estimated  to 
be  r"o=  5.4  km,  which  is  within  10%  of  the  true  source  range. 

The  approach  used  for  localization  of  a  single  source  can  also  be  used  for  localization  of  multiple 
uncorrelated  random  noise  sources  in  an  ocean  waveguide  in  both  range  and  bearing.  Here  we 
consider  the  case  where  there  are  3  uncorrelated  noise  sources  SI,  S2,  S3,  as  shown  in  Fig.  4.  The 
source  SI  is  at  range  rl  =  8.2  km  and  bearing  $\theta_l  =  45A\circ$.  The  other  two  sources  S2  and  S3 
are  at  the  same  bearing  $\theta_2  =  \theta_3  =  60A\circ$,  and  their  ranges  are  r2  =  3  km  and  r3  =  10 
km,  respectively.  We  assume  that  these  three  sources  have  the  same  source  power  spectrum  of  0  dB  re 
1  pPa2  /  Hz  in  the  390  to  440  Hz  frequency  band,  but  they  are  uncorrelated  with  each  other. 
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Figure  4:  Top  view  of  the  geometry  for  the  multiple  source  localization  example.  All  3  sources  are 
at  30-m  depth,  and  the  receiver  array  depth  is  50  m.  The  receiver  array  aperture  is  150  m. 

The  cross-correlated  intensity  field  IBo  (s,x  )  for  this  multiple  source  scenario  is  shown  in  Fig.  5.  The 
cross-correlated  intensity  in  Fig.  5  exhibits  two  distinct  intensity  fields  when  x  <  0,  one  at  the  bearing 
of  SI,  and  the  other  at  the  bearing  of  S2  and  S3.  When  x  >  0,  the  field  for  S2  and  S3  splits  into  two 
distinct  fields,  since  their  ranges  from  the  receiver  array  are  different.  The  cross-correlated  intensity 
field  for  SI  when  x  >  0  shows  the  same  pattern  as  that  of  the  single  source  scenario. 


Cross-correlation 


Figure  5:  The  cross-correlated  intensity  field  IBo  (s,t  )  for  the  source-receiver  geometry  shown  in 
Fig.  4.  The  two  black  solid  lines  overlain  in  the  t<0  domain  marks  sinO  and  sinO ,  respectively.  The 
three  black  solid  lines  in  the  12  t>  0  domain  are  the  beamformer  migration  lines  for  sources  SI,  S2, 

and  S3  in  terms  of  the  reduced  travel  time. 

We  now  can  apply  the  image  transform  method  previously  introduced  for  localization  of  single  source, 
and  localize  all  three  sources  in  range  and  bearing  without  ambiguity.  The  detailed  procedures  and 
estimation  results  are  provided  in  Ref.  [7]. 
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IMPACT/APPLICATIONS 


The  array  invariant  method  developed  from  Fiscal  Year  2006  effort  has  significant  advantages  over 
existing  source  localization  methods  for  practical  source  localization  scenarios.  This  is  because  the 
array  invariant  method  does  not  require  a  priori  knowledge  of  the  environmental  parameters,  nor  does 
it  require  extensive  computations. 

The  array  invariant  method  has  tremendous  potential  application  in  shallow-water  surveillance 
missions  and  anti-submarine  warfare.  The  method  enables  instantaneous  and  simultaneous  localization 
of  multiple  broadband  noise  sources  in  a  littoral  environment  using  towed  arrays  from  surface  ships  or 
submarines. 
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